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We report measurements of the polarization of W bosons from top-quark decays using 2.7 fb~^ 
of pp collisions collected by the CDF II detector. Assuming a top-quark mass of 175 GeV/c^ three 
measurements are performed. A simultaneous measurement of the fraction of longitudinal (/o) and 
right-handed (/+) W bosons yields the model-independent results fo = 0.88±0.11(stat) ±0.06(syst) 
and /+ — —0.15 ± 0.07(stat) ± 0.06(syst) with a correlation coefficient of -0.59. A measurement of 
fo (/+) constraining /+ (/o) to its standard model value of 0.0 (0.7) yields fo = 0.70 ± 0.07(stat) ± 
0.04(syst) (/+ — —0.01 ± 0.02(stat) ± 0.05(syst)). All these results are consistent with standard 
model expectations. 

PACS numbers: 14.65.Ha, 14.80.Cp, 13.85.Ni, 13.85.Qk, 14.70.Dj 
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The top quark is the most massive fundamental par- 
ticle observed by experiment [ij. Due to its large mass, 
in the standard model (SM) the top quark decays be- 
fore forming a bound state via the charged current weak 
interaction into a boson and a b quark Q, with a 
branching fraction above 99% Q • This provides a unique 
opportunity to study the properties of a "bare" quark. 
In particular, the V — A structure of the weak interac- 
tion can be tested by reconstructing the polarization of 
the W'^ boson from top-quark decay. In the standard 
model (SM) at tree level [J], the boson is expected 
to have longitudinal polarization /o = 0.703, left-handed 
polarization /_ = 0.297, and right-handed polarization 
/+ = 3.6x 10^** for a top-quark mass rrit = 175 GeV/c^, a 
ly-boson mass Mw = 80.413 GeV/c2 @, and a 6-quark 
mass mi, = 4.79 GeV/c^ In the limit of irib — 0, 
/o = TOj/(2m^ + TOj) and /+ = 0. The higher-order 
QCD and electroweak radiative corrections modify these 
predictions at the 1-2% (relative) level Q. In beyond- 
the-SM scenarios, significant deviations from the SM ex- 
pectation are possible due to the presence of anomalous 
couplings Q in the tWb vertex. Measurements of W- 
boson polarization and single top-quark production to- 
gether set constraints on the anomalous coupling vector 
and tensor form factors Q- 

In this Letter, we measure the polarization of the W 
boson from top-quark decay. We assume the ti pro- 
duction mechanism is in agreement with the SM, and 
we study a data sample enriched in tt W^bW^b -> 
tvbqq'h events where one of the W bosons decays hadron- 
ically and the other leptonically. We apply a likeli- 
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hood technique based on the theoretical matrix elements 
for both the dominant signal process, qq — ^ it, and 
the main background process, inclusive production of 
W-|-jets. This technique was first developed for the mea- 
surements of top-quark mass and /o, constraining /+ to 
its SM value [8|, and utilizes the kinematic and topo- 
logical information from the event through integrations 
over poorly known parton-level quantities. We express 
the matrix element in terms of the VF-boson polariza- 
tion fractions and the cosine of the angle 9* between the 
momentum of the charged lepton or down-type quark 
in the VF-boson rest frame and the momentum of the 
W boson in the top-quark rest frame. Therefore we 
extract information on the VF-boson polarization from 
both the leptonic and hadronic V(^-boson decays. Pre- 
vious CDF measur ements[i, ^ used only information 
from the leptonic decay via the transverse momentum 
ll| of the lepton, cos 6* for the leptonically decaying 



Pt 

W boson, or the invariant mass of each lepton and 6-jet 
pair, which is a function of cos 9* . While the information 
from the hadronic W-boson decay carries a sign ambigu- 
ity in cos 9* since we are unable to identify the down-type 
quark jet its inclusion still improves the sensitivity to the 
/o polarization fraction. The analysis described in this 
letter improves the statistical sensitivity on /o by 20% 
relative to the best previous CDF measurement Q for 
the same event sample. The latest DO measurement also 
utilizes information from both the leptonic and hadronic 
IF-boson decays using a likelihood fit to the cos 6'* dis- 
tributions jTij . 

We report measurements of the Ty-boson polarization 
for three different hypotheses of top-quark decay: (1) 
model-independent with simultaneous measurement of 
/o and /+; (2) anomalous tensor couplings with mea- 
surement of /o for fixed /+=0; and (3) anomalous right- 
handed couplings with measurement of for fixed /o = 
0.70. 

The polarization fractions are determined by maximiz- 
ing the likelihood L with respect to /o, /+, and the frac- 
tion of events consistent with the tt signal hypothesis, 



N 



L(/o,/+,a) = nc« 



+ (1 - C,) 



Pbjx) 
(Mx)) 



Here N is the number of observed events, and (As) and 
(Ab) refer to the average acceptances for ti and W-|-jets 
background events, respectively. The dependence of the 
signal acceptance on the polarization fractions is included 
in {As). The signal probability Pg and b ackg round prob- 
ability Pb densities are constructed as in [13[ by integrat- 
ing over the appropriate parton-level differential cross 
section, da{y)/dy, convolved with proton parton distri- 
bution functions (/(^i) and f{q2)) and detector resolu- 
tion effects relating a set of observed variables x to cor- 
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responding parton- level quantities y, W{x^y): 



d<j{y) 

dy 



f{qi)fiq2)dqidq2W(x, y)dy. 



Here partons are identified with the four highest trans- 
verse energy Et 11| jets, and all the corresponding jet- 
parton permutations are considered. 

The signal differential cross section uses the leading- 
order matrix element of the qq — > ti process 14|, ex- 
pressed in terms of cos 6* and polarization fractions: 



WY = ^Fi.pFhadi'i - /?'sin^0,t), 

where gs is the strong coupling constant and Oqt describes 
the angle between the incoming parton and the top quark 
in the rest frame of the incoming partons, and /3 is the 
speed of the top quarks in the same rest frame. The 
factors Fiep and Fhad correspond to the top quarks with 
a leptonic and a hadronic ly-boson decay, such that: 



F, 



lep 



27r 5h'™L.„p*2 , op*,^ 



ml) 



3 mtTt 

(/+.^(i + cosr)2-H/o-^(i-cos2r) 
+(i-/o-/+)-|(i-cosr)2). 



Here gw is the weak coupling constant, mjf^ is the invari- 
ant mass of the lepton and neutrino, Tt is the width of 
the top quark, mt and mb are the masses of the top quark 

and b quark, respectively, and = ' — ^. The 

hadronic factor F^ad is similar, with the exception that 
we do not distinguish between up-type and down-type 
quarks from VF-boson decay and use the average Fhad 
related to the two permutations. The background differ- 
ential cross section uses the sum of matrix elements for 



W-t-jets from the VECBOS [151 Monte Carlo (MC) gener- 
ator. 

The measurement is based on a data set with an inte- 
grated luminosity of approximately 2.7 fb~^ acquired by 
the Collider Detector at Fermilab (CDF H) [3l from pp 
collisions at ^/s = 1.96 TeV. The data used are collected 
using high-pT central (pseudorapidity [llj |?7| < 1.1) 
electron and muon triggers, a high-p^ forward (1.2 < 
\ri\ < 2.0) electron trigger, and a trigger that requires 
large missing transverse energy [ill with either an 
energetic electromagnetic cluster or two separated jets 
(^T+jets) [l3|. The ^T+jets trigger is used to select 
additional events with high-py muons, which are not se- 
lected by the lepton triggers. 

Candidate events for the lepton plus jets final state 



{ti — W~^bW~b — )• ivhqq'h) are selected to have a single, 
isolated electron or muon candidate with Et > 20 GeV, 
large imbalance in transverse momentum in the event 
{^T > 20 GeV) as expected from the undetectable neu- 
trino, and at least four jets with Et > 20 GeV. Jets 
are reconstructed using a cone algorithm with radius 
AR = 0.4 in 77 (j) space, and their energies are cor- 
rected for non-uniformities in the calorimeter response 
as a function of jet 77, multiple pp interactions, and the 
hadronic jet energy scale of the calorimeter [l^l . Of these 
jets, we require at least one to have originated from a 
b quark by using an algorithm that identifies a long- 
lived b hadron through the presence of a displaced ver- 
tex (b tag) [3]. Backgrounds to the tt signal arise from 
multi-jet QCD production (QCD), VF-boson production 
in association with jets (W^-l-jets), and electroweak back- 
grounds (EWK) composed of diboson {WW, WZ, ZZ) 
and single top-quark production. The VF-l-jets back- 
ground includes 6-flavor jets (M^-|-hf) as well as light fla- 
vor jets incorrectly identified as b jets (M^-flf). 

A detailed description of the background estimation 
can be found in Ref. [2^ . Table U shows the expected 
sample composition assuming a tt cross section of 6.7 pb. 
There are overlapping events between those collected by 
high-pT lepton triggers and ^r+jets trigger. These over- 
lapping events are accounted in the central e/^ and for- 
ward e categories, and are explicitly vetoed from ^T+jets 
category. 



TABLE I: Number of expected and observed events in 2.7 
fb^^ assuming a ti cross section of 6.7 pb. 



Process 


Central 


Forward 


^T+jctS 




e, ^ 


e 




tt 


478 ± 66 


58 ± 8 


134 ± 19 


W+M 


71 ± 22 


13 ± 9 


19 ± 6 


W^-l-lf 


23 ± 6 


5 ± 7 


6 ± 2 


EWK 


17 ± 10 


3 ± 1 


5 ± 3 


QCD 


28 ± 22 


46 ± 37 


1 ± 1 


Total expected 


616 ± 74 


125 ± 40 


165 ± 20 


Observed 


650 


136 


178 



The HERWIG [21| MC generator is used to model the ti 
signal events with mt = 175 GeV/c^. For estimation of 
various systematic uncertainties and background model- 
ing MC samples are created using the pythia 2^ gener- 



24l | with PYTHIA or 



ator, and ALPGEN [23[ or madevent 
HERWIG supplying the parton shower and fragmentation. 
The QCD background is modeled using data control sam- 
ples. The signal and background modeling has been ex- 
tensively checked. Figure [1] compares the observed data 
and the MC-predicted distributions of different kinematic 
variables. We have validated the background model by 
studying a high-statistics control sample of VF-|-jets can- 
didates extracted by vetoing events containing 6-tagged 
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jets. 
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FIG. 1: Comparison of four kinematic variables for data 
and simulation for different W polarization fractions: solid, 
dashed and dotted histograms correspond to (/o, /+) values 
of (0.7, 0.0), (0.88,0.0) and (0.7, 0.3), respectively. Plotted 
are (a) leading jet pr, (b) lepton pt; and for the reconstruc- 
tion chosen as most likely by the per-event likelihood (c) the 
invariant mass of the pair of light quark jets from the hadron- 
ically decaying W boson and (d) the cos 9* of the leptonically 
decaying W boson. 



cesses. Even though likelihood can be calculated only 
for the physical values of fo and /+ , after calibration the 
corrected measured values can be slightly outside their 
physical ranges. 

The robustness of the fitting procedure over all physi- 
cal values of (/o,/+) has been tested with simulated ex- 
periments, using the number of observed data events and 
the sample composition of table HI In all cases, the fit is 
unbiased. Near the physical boundaries, we find that the 
statistical uncertainty is underestimated by as much as 
a factor of 1.5. We apply a correction to the statisti- 
cal uncertainty in these regions. Assuming the SM, the 
expected statistical uncertainties after all corrections for 
the simultaneous measurement are ±0.116 and ±0.074 
for /o and /+, respectively. 

Various sources of systematic uncertainty affecting the 
measurement arc summarized in tabic |lTl The leading 
sources of systematic uncertainty arise from MC model- 
ing of initial and final state radiation (ISR/FSR), choice 
of parton distribution functions (PDFs), choice of par- 
ton shower model, uncertainties on the measured jet en- 
ergy, and the background shape and normalization. The 
method-related uncertainty includes propagating the un- 
certainty on the fit parameters of the response curves, 
including their correlations. All systematic uncertain- 
ties arc determined by performing simulated experiments 
in which the systematic parameter in question is varied, 
the default method and calibrations are applied, and the 
shifts in the mean measured polarization fractions are 
used to quantify the uncertainty. All shifts are evaluated 
at the SM helicity fraction. 



We calibrate the results of the likelihood fit using 
the simulated tt and background samples, and the sam- 
ple composition of table ID For the simultaneous mea- 



surement of fo and /- 



+ 1 



find our estimate /o,i 



IS 

related to the true value of fo by /o,m = (0.88 ± 
0.02)/o + (-0.12 ± 0.01) and our estimate of /+,,„ is 
related to the true value of /+ and fo by /+.m = 
(1.26 ± 0.01)/+ + (0.17 ± 0.02j/o + (0.06 ± 0.01)^ We 
use these two equations and the measured polarization 
fractions to extract the true polarization fractions. For 
our measurement of fo with /+ = 0, we find our esti- 
mate /o,m = (1.15±0.04)/o + (-0.09±0.02),andfor our 
measurement of /+ with fo = 0.7, we find our estimate 
/+a„ = (1.17 ± 0.05)/+ + (0.01 ± 0.01). The uncertain- 
ties on the coefficients of the calibration functions are 
accounted as the method-related systematic uncertain- 
ties, which cover possible biases due to the calibration 
procedure. The differences between /o,m or /+.m values 
and the corresponding true values are due to the fact that 
the signal and background probabilities in the the likeli- 
hood do not model events with extra jets from initial and 
final state radiation (ISR/FSR), or incorrect jet-parton 
assignment where a ISR/FSR jet is selected as one of 
the leading four jets, or all the different background pro- 



TABLE II: Summary of systematic uncertainties. 



Source 


A/o 


A/+ 


A/o 


A/+ 








simultaneous 


ISR/FSR 


0.020 


0.018 


0.020 


0.021 


PDF 


0.024 


0.013 


0.009 


0.016 


JES 


0.018 


0.017 


0.004 


0.012 


Parton shower 


0.012 


0.008 


0.031 


0.017 


Background 


0.009 


0.038 


0.042 


0.039 


Method-related 


0.010 


0.005 


0.024 


0.024 


b-tag SF 


0.004 


0.002 


0.002 


0.002 


Total 


0.041 


0.048 


0.062 


0.057 



For the simultaneous measurement of fo and /+, we 
exclude the events from the forward electron trigger as 
this significantly reduces the systematic uncertainty from 
the background model. With 828 events and after all 
corrections, we measure 

fo = 0.879 ±0.106(stat)± 0.062(syst) 
/+ = -0.151 ±0.067(stat)± 0.057(syst). 
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The statistical correlation between /o and /_)_ is p = 
-0.59. We estimate a shift of t(0.010 ± 0.005) in /o 
and ±(0.017 ± 0.003) in /+ per ±1 GeV/c^ shift in the 
top quark mass from the central value of 175 GeV/c^. As 
the central value is unphysical we have elected to ensure 
coverage by applying the Feldman-Cousins method 
to obtain the confidence level (C.L.) intervals shown in 
Figure EI 

Fixing /+ = and with 964 events, we measure af- 
ter all corrections /o = 0.701 ± 0.069(stat) ± 0.041(syst). 
Fixing /o = 0.70, we measure after all corrections /+ = 
-0.010 ± 0.019(stat) ±0.049(syst) and find /+ < 0.12 at 
95% C.L.. We estimate a shift of ±(0.011 ± 0.003) in /o 
and ±(0.013 ± 0.002) in /+ per ±1 GeV/c^ shift in the 
top-quark mass from the central value of 175 GeV/c^. 



0.5 









— 68% C.L. 




90% C.L. 







0.5 1 

fo 

FIG. 2: Contours in the (/o,/+) plane indicating 68% 
and 90% C.L. intervals determined using Feldman-Cousins 
method. Note that the coverage is correct although the cen- 
ter of the contours is not at the measured value obtained after 
calibration. 

In summary, we have measured the polarization of the 
W boson in top-quark decays using a matrix-element 
method in 2.7 fb~^ of CDF II data. Our results are con- 
sistent with the SM. This result improves the combined 
statistical and systematic precision on both the model- 
independent and model-dependent determinations of the 
longitudinal polarization /p by a factor of 1.3 compared 
to the best previous measurement Q for a 1.4 times in- 
crease in luminosity. 



Wc thank the Fcrmilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda- 
tion; the A. P. Sloan Foundation; the Bundesministerium 
fiir Bildung und Forschung, Germany; the World Class 
University Program, the National Research Foundation 
of Korea; the Science and Technology Facilities Coun- 
cil and the Royal Society, UK; the Institut National de 
Physique Nucleaire et Physique des Particules/CNRS; 
the Russian Foundation for Basic Research; the Minis- 
terio de Ciencia e Innovacion, and Programa Consolider- 
Ingenio 2010, Spain; the Slovak R&D Agency; and the 
Academy of Finland. 



[1] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 74, 
2626 (1995); S. Abachi et al. (DO Collaboration), Phys. 
Rev. Lett. 74, 2632 (1995). 

[2] Charge-conjugate modes are included implicitly through- 
out this paper. 

[3] C. Amsler et al. (Particle Data Group), Phys Lett. B667, 

1 (2008) and 2009 partial update for the 2010 edition. 
[4] J. Aguilar-Saavedra et al, Eur. Phys. J. C 50, 519 (2007). 
[5] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 

99, 151801 (2007). 

[6] M. Fischer et al, Phys. Rev. D 63, 031501 (R) (2001); H. 
S. Do et al, Phys. Rev. D 67, 091501(R) (2003). 

[7] V. M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
102, 092002 (2009). 

[8] V. M. Abazov et al (DO Collaboration), Nature 429, 
02589 (2004); V. M. Abazov et al (DO Collaboration), 
Phys. Lett. B617, 1 (2005). 

[9] T. Aaltonen et al (CDF Collaboration), Phys. Lett. 
B674, 160 (2009). 
[10] A. Abulencia et al. (CDF Collaboration), Phys. Rev. 
Lett. 98, 072001 (2007); A. Abulencia et al (CDF Col- 
laboration), Phys. Rev. D 75, 052001 (2007); A. Abulen- 
cia et al (CDF Collaboration), Phys. Rev. D 73, 111103 
(2006); D. Acosta et al (CDF Collaboration), Phys. Rev. 
D 71, 031101 (2005). 
[11] We use a cylindrical coordinate system where the z axis 
is along the proton beam direction and 6 is the polar 
angle. Pseudorapidity is -q = — lntan(0/2), while trans- 
verse momentum is pr = \p\ sinO, and transverse energy 
is Et = EsinO. Missing transverse energy, ^t, is defined 
as the magnitude of — J]]^ E^fii, where fii is the unit vec- 
tor in the azimuthal plane that points from the beam line 
to the i*'' calorimeter tower. 
[12] V. M. Abazov et al (DO Collaboration), Phys. Rev. Lett. 

100, 062004 (2008). 

[13] A. Abulencia et al (CDF Collaboration), Phys. Rev. 

Lett. 99, 182002 (2007). 
[14] G. Mahlon and S. Parke, Phys. Lett. B411, 173 (1997); 



8 



G. Mahlon and S. Parke, Phys. Rev. D 53, 4886 (1996). 

[15] F. A. Berends, W. T. Giele, and H. Kuijf, Nucl. Phys. 
B321, 39 (1989). 

[16] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
031101(R) (2005). 

[17] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 
74, 072006 (2006); T. Aaltonen et al. (CDF Collabora- 
tion), Phys. Rev. Lett. 100, 211801 (2008). 

[18] A. Bhatti et al., Nucl. Instrum. Methods A 566, 2 (2006). 

[19] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
052003 (2005). 



[20] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 

97, 082004 (2006). 
[21] G. Corcella et al, J. High Energy Phys. 01, 10 (2001). 
[22] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 

(2001). 

[23] M. Mangano et al., J. High Energy Phys. 07, 001 (2003). 
[24] J. Alwall et al., J. High Energy Phys. 09, 28 (2007). 
[25] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 
(1998). 



